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Abstract 

This work considers various families of quantum control landscapes (i.e. objective functions for 
optimal control) for obtaining target unitary transformations as the general solution of the controlled 
Schrodinger equation. We examine the critical topology of the kinematic landscapes Jf(U) = ||(t/ — 
W)Af and Jp{U) = II^H" - | TT{AA^W^U)f defined on the unitary group U(iV). The parameter matrix 
A is allowed to be completely arbitrary, yielding an objective function that measures the difference in 
the actions of U and the target on a subspace of state space. The analysis of this function includes 
a description of the structure of the critical sets of these kinematic landscapes and characterization of 
the critical points as maxima, minima, and saddles. In addition, we consider the question of whether 
these landscapes are Morse-Bott functions on U(A^). These results are then used to deduce properties 
of the critical set of the corresponding dynamical landscapes. Finally, landscapes based on the intrinsic 
(geodesic) distance on U(A^) and the projective unitary group PU(A'') are considered. 

1 Introduction 

An important application of quantum optimal control theory is the generation of target quantum logic gates 
for quantum information processing. The goal of such optimal control is to arrange the dynamics such that 
the desired logical gate is realized as the final time unitary evolution operator, which is the general solution 
of the controlled Schrodinger equation. Additionally, there may not be a single target gate, but a family 
of them. For example, since the global phase is not observable, the goal may be any unitary operator that 
is equivalent to the target gate up to global phase. Likewise, in some cases only a subspace of the Hilbert 
space of states may be used for the quantum register, so that any unitary propagator should be acceptable 
that acts as the target gate on that subspace. In contrast to other quantum control problems, for example 
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the maximization of a quantum mechanical observable, there is no unique or natural choice for the objective 
function against which the optimization is performed. Indeed, any smooth function J : V{N) — >■ M with a 
global minimum at the target unitary gate or gates is a candidate objective for the unitary problem. But, 
as we will see, some choices may have more favorable convergence and other properties. 

Let K denote the space of admissible control functions. For the present analysis, K will be L'^(M+;M), 
where T is some fixed final time over which the controlled dynamics take place. Let C^^^ possess the real 
Hilbert-Schmidt inner product {A, B) = ^Ty{A^B), and let U(iV) C C^""^ denote the x iV unitary group 
endowed with the Riemannian metric induced by this inner product. Also let Vr : K — > U(Af) denote the 
map, defined implicitly by the Schrodinger equation in the dipole approximation 

ih^{tM) = {HQ- ti£{t))U{tM) U{to,to)^l (1) 

such that Vt{£) = U{T,0)[£] is the unitary propagator at time T for the control field £. Finally, for any 
candidate objective function J : U(A^) — > R (the "kinematic landscape"), let J : K ^ M (the "dynamical 
landscape") be the composition J = J oVr- Then 

gradJ(£) - (d£W)*(grad J(1/t(5))), (2) 

where (d^V^)* is the operator adjoint of the differential dgVr- Much of the important information about the 
nature of the gradient flow of J is embodied in the critical points of this landscape, i.e. those flelds £ e K 
for which grad J(£) = 0. Any £ G K such that dgYr is full rank and grad J(Vt(i^^)) = (so-called "regular" 
critical points) will satisfy the condition. There may be other critical points where dfVr is rank-deficient 
and grad j(VT(f)) may or may not be zero ("singular" points). Consideration of such singular points is 
important for a complete understanding on the dynamical control landscape. However, since Vr is a highly 
nonlinear map from an infinite dimensional space to a finite dimensional space, singular points are expected 
to be rare and will not be considered in the present analysis. Singular points and their role in quantum 
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control landscapes have only recently begun to be analyzed 

The critical points of the kinematic landscape having been identified, they may then be characterized as 
local maxima, local minima, and saddles. As has been demonstrated for other classes of kinematic quantum 
control landscapes 0, [ll|, [l2, 16-19, 22, 23 1, the landscapes considered in this work will turn out to have 
global maxima and minima, but no other local extrema. Moreover, the critical sets will be shown generally 
to comprise disjoint submanifolds, and these submanifolds are nondegenerate in the Morse-Bott sense. In 
other words, the null space of the Hessian of J and the tangent space of the critical submanifold coincide at 
each critical point U G V{N). This condition identifies the kinematic landscape as a Morse-Bott function, 
which is interesting for at least two reasons. First, certain results about the convergence of the gradient 
fiow may be proved for Morse-Bott functions, in particular that (on a compact manifold) the gradient flow 
always converges to a critical point [7|. Second, the identification of the null space of the Hessian and the 
tangent space of the critical submanifold is important for certain numerical methods, such as second order 
D-MORPH ^], that are designed to explore the critical sets. 
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Several classes of landscapes for generating target unitary transformations will be considered. They 
include Jf{U) = \\{U - W)Af and Jp{U) = - \Ti {AA'^W^U)]^ for some fixed target W e U(iV) 
and some arbitrary fixed A G C^^^, as well as the corresponding landscapes using the intrinsic (geodesic) 
distance on the unitary group V{N) and the projective unitary group PU(iV), rather than the norm (Eu- 
clidean) distance in Jp and Jp. Landscapes of the form Jp and Jp have been studied in the past 
The present paper extends these various works by broadening the families of landscapes under considera- 
tion, describing the structure of the critical submanifolds, and directly addressing the issue of Morse-Bott 
nondegeneracy of the critical submanifolds. 

The paper is organized as follows. Sections[5]and[3]describe the critical points of the kinematic landscapes 
Jp and Jp. In Section |4l these results are related back to the dynamical landscapes. The two additional 
landscapes based on geodesic distance are presented and analyzed in Section [S] The overall results are 
summarized in Section [51 Two appendices are included which provide a proof of the infinite Frechet dif- 
ferentiability of the control-to-propagator map Vp and a derivation of the gradient of the geodesic distance 
landscape on V{N). 



2 Kinematic Critical Point Analysis of Jf{U) = \\{U - W)A\\'^ 

For now we put aside the dynamical component of the map and focus just on the critical point analysis of the 
kinematic map Jp{U) ^ \\{U ~ W)A\\'^ = 2||A|p - 2^Tr{AA'<W''U). In contrast to the landscape Jp that 
wiU be the subject of the next section, the value of Jp{U) depends upon the global phase of U. However, if 
AA^ has a null space, then Jf{U) is invariant to the action of U on that null space. This freedom allows for 
the possibility that A is a projection matrix as discussed in [15|. 



2.1 Critical Point Identification 

Let D e U(Af) be a unitary matrix that diagonalizes AA'^ such that D^AA^D = with Qjj = ujj and 
ujf > a;| > • • • > . Let uif > ■ ■ ■ > uj^ > ojq = denote the distinct eigenvalues of AA^ , and let 
{ni, . . . , n^, hq} be the degeneracies of these eigenvalues. While ni, . . .rii^ are all strictly positive integers, 
riQ may be either positive or zero. This distinction is made due to the special significance of the zero 
eigenvalues of AA^ in the analysis that follows. Suppose that D E V{N) is any unitary matrix such that 
Dn^D^ = AA^ = Dn^D^. Then {D^ b)^'^ {D^ 1))^ = n^, so that D lies in the stabilizer subgroup of n^, 
i.e. D^iJ) e g := U(ni) ® • • • ® U(n«,) ® U(no) C V{N) and D e DQ. In other words, DQ is the set of all 
unitary matrices {D} such that D^AA^D = fl^. 

Now, the differential of Jp, djjJp ■ T;7U(iV) — > M, is given by 

duJp{SU) = -2mr{AA^W^dU) = {U AA^ U - W AA\ SU) , (3) 

where the last step included a projection of ~2WAA^ into the tangent space TijlJ{N). Therefore, 

gra.d Jp{U) = UAA'fW^U - WAA^ (4) 
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and a critical point of JpisnU such that UAA^W^U = WAA^. Let X = W^U £ U(iV), and let X = D^XD. 
Then the critical point condition becomes X^Q^ = V,^X. This result also implies that ATi^ — Vl'^X . Hence 
for any i and j, if ^ then Xij = 0, and if = w| 7^ then A^^j = So, X is unitary and block- 
diagonal with block sizes {rii, . . . , n^, no} and the first k blocks are Hermitian. In particular, this implies that 
[X, AA"^] = 0, so that X and AA"^ may be simultaneously diagonalized by some DT G DQ. Every critical X 
can then be written as A" = D{tkt^ ®Xq)D\ where f e G U(ni)® • • •©U(n„) C U(7V), A = Ai® . . . A^, 
Aj = — Iiy. ffi Irii-i/i G U(ni), < i^i < Tii, and G U(no). Only the number Vi of negative I's and Ui — Vi of 
positive I's is important in defining A^, not their order, since F can apply any permutation to the elements 
of Ai. The order above is chosen for convenience. Moreover, every such X ~ Z)(rAr^ Ao)-D^ with F G ^, 
A = Ini-i/i • • • ® ^Ii^K ffi In^-i/^, and Ag G U(no) corresponds to a critical point of Jp with critical 
value Jf{U) ~ 4:^21=1 ^i^l- Thus, the set of all critical points is 

Crit( Jf) = {u = WD{tKt^ ® Xo)D^ : 

f G ^, A = -I^, e Ini-v^ • • • © © In^-i'., and ib G U(no)|. (5) 

Let $ : (? X U(7V) U(iV) be the Lie group action $(F, V) = TV{t^ © !„(,), where f is the projection 
of F into Q. For a fixed V G U(iV), the set of all matrices reached by the action is the orhit through V: 
Oih^iV) ^ {$(F, y) : F G G}. Similarly, for a fixed V G U(7V), the set of aU F's in G that keep V fixed is 
called the stabilizer (or isotropy) subgroup of G'- Stab$(y) = {V £ G ■ $(F, = V}. Since is a compact 
Lie group, Orb$(y) is an embedded submanifold of U(A^) for each V G U(7V). Furthermore, Orb$(y) 
is diffeomorphic to the connected diffcrcntiable manifold G / Stah^{V). In particular, dim(Orb$(y)) = 
dim(^) — dim (Stab$(y)). Details on this orbit-stabilizer result may be found, for example, in [ist or 
Ch. II, Ex. C.5]. 

For £L given A — — li/-^ © 1^^^—^^^ • • • — Ii^^ © I^i^— i/^, let us now characterize Orb$(A -lino)- Observe that 
Stab$(A©I„o) = U(i/i)©U(ni-j/i)©- • •©U(:/«;)©U(n„-j/„)©I„o. The dimension of ^ is dim(t;) ^J2i=o'^h 
while the dimension of Stab$(A © I„o) is dim ( Stab$(A © Ino)) = J2i=i ^1 + ("-i ~ ^0^- Hence, we obtain 

nu/Affif ^ U(ni)©---©U(nJ©U(no) 

© i„o J - ^ - ^i) © • • • © U(^.) © U(n. - © I„„ ^^^^ 

~ Gr,, (C"i) © • • • © Gr,^ (C"») © U(no) (6b) 

dim ( OrbcE,(A © I„ J) ="0 + 51"?" ~ ~ = "o + 2 XI ^^("« ~ '^')' (^c) 

where Gr^(C") denotes the complex Grassmannian manifold of j/-dimensional complex subspaces of C". 
Therefore if AA^ has k distinct non-zero eigenvalues with degeneracies {ni, . . . ,71,^}, then Jp has exactly 
M = nr=i('^« + ^) connected critical submanifolds. Note that in the case of complete degeneracy of AA''' 
(e.g. A — I), K = 1 and ni = N, implying that there are N + 1 critical submanifolds with dimensions 
0, 2{N — 1), 4{N — 2), 6(iV — 3), . . . for I'l = 0, 1, 2, 3, ... . In the case of complete non-degeneracy and non- 
singularity of AA^ , K = N and ni = 1 for all i, so there are 2^ critical submanifolds, all of dimension 0, i.e. 
isolated points. 
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2.2 Hessian Analysis 



Turning to the question of the signatures of these critical points, we extend the gradient vector field grad Jp 
in the obvious way to all of C^^^ and differentiate to find du grad Ji.((5C/) = 5U AA'^W'^U + UAA^WUU. 
Projecting this onto the tangent bundle of U(-/V) gives 



Hessjp,t/((5J7) := V su grad Jf(C/) 

= ]^{5UAA^W^U ~ WAA^SU^U + UAA^W^SU - USU'^WAA'') 
= -{SUAA^W^U + WAA^U^SU + UAA^W^SU + SU U^WAA^), 



(7a) 
(7b) 
(7c) 



using the fact that any tangent vector SU € TijV{N) satisfies SW = —WSUW, and where Vsu denotes 
the covariant derivative in the direction SU [3]. At a critical point, grad Ji?(J7) = 0, so that W'^U and WW 
both commute with AA'' , and also W^'UAA'' = AA^WW = Dn^{tAt'' ® Xo)D^ = DQ^{tAr^ ® 0)£)t = 
WWAA'f = AAWW. Then the Hessian becomes 



Hess 



Jf,U 



{6U) = SUAA^W^U + UAA^W'fSU. 



(8) 



Suppose that U G V{N) is a critical point of Jp, i-e. U = WD(tAt'^ ® Xn„)D^ for some f e ^ and 
A = -I^, •••®-I^^®I„^_i,^ and i^o G U(no). Furthermore, let F := (ft ®I„j£)t[/tj[/£)(f®I„J 

so that Y is skew-Hermitian. Then Ress j^{SU) = jSU if and only if Yn'^{A ® I„J + (A ® I„o)^ = 7^- 
For the null space of Hessj^, we want to solve this equation for 7 = 0. This equation is solved in u(A^) if 
and only if Y is of the form 



Y = 



Qi 
-Ql 



Qk 



Qo 



(9) 



where Qi G for i = 1, . . . , k and Qo G u(no) are arbitrary. These F's span a dimension d 



- 2 ^ t'i(rii — Vi) (real) subspace of u(A^). Hence the nullity of Hessj^ is Ao 



-2J2'^i{ni- Vi)- Since 



this is the same as the dimension of the corresponding critical submanifold, and since the tangent space of 
the critical submanifold must lie in the null space of the Hessian, these two subspaces of TuV{N) must be 
identical. Therefore Jp is Morse-Bott for all A matrices. 

To conclude this kinematic analysis, we can compute the numbers of positive and negative eigenvalues (i.e. 
Af+ and Af- , respectively) of Hessj^^ at a critical point U. Let Ai = A^^ for i = 1, . . . , N — uq and Ai = 1 for 
i = N-na + l,...,N. Then rfi2(A®I„J + f^^(A®I„jy = jY can be expressed as {ujfXj+ujfXi-j)Yij = 
for 1 < i, j < N . Then 7 is an eigenvalue of Hessj^ if and only if 7 = i^j Aj + ujfXi for some i and j. Now, 
w^Aj + ujf Xi < whenever: 
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• Ai = Xj = —1, which happens for (X]r=i ^0 ordered pairs 

• \i = —1, Aj = +1 and uf > uj'j, which happens for J2k>i>i ^k{ni — vi) + nQ '}2d=i ordered pairs (i, j), 
or 

• A, = +1, Xj = -1 and io'j < uj'^, which happens for J2i<k<ii'^k - Vk)vi + rio Ya=i = J2k>i>i ^^{ni - 
vi) + riQ X]r=i ordered pairs (i, j). 

Since each ordered pair (i, j) contributes one (real) degree of freedom in K, there are A/"- = (X]r=i '^O^ + 
2Sfc>i>i ^fc('T-i — i^i) + 2no X]r=i negative eigenvalues of Hessj^. Likewise, Aj + tof Xi > whenever 

• Ai = A, = +1 and at least one of ujf and is nonzero, which happens for 
(iV - no - ^i) + 2"-o (A^ - no - Z]fe=i ^k) ordered pairs (i, j), 

• A; = —1, Xj = +1 and ujf < which happens for X]i<fe<i ^k{ni — vi) ordered pairs («, j), or 

• Xi = +1, Aj = -1 and ujf > a;|, which happens for Y,k>i>i(''^k ^ Vk)vi = X]i<fc<; '^*:('^' ^ ^i) ordered 
pairs (i, j). 

Therefore, the positive eigenvalues number 

N+^iN-no-Y^vA +2no (iV-no-^i/fc J +2 ^ z.fe(nj-z/;) (10a) 

V 1=1 / \ k=l / l<k<l 

= 7V2-n2+ (^z/,)'-2iV^j/, + 2 ^ Vk{ni-vi). (10b) 

/=i /=i i<fc<i 

It is easy to see that M- + N+ + No = N'^ as expected. Furthermore, we find that M+ = if and only if 
Vi — rii for alH = 1, . . . , K, i.e. only at the global maximum {U — — M^£)(I©<Yo)D^ : Xq e U(no)}. Likewise 
M- = if and only if = for all « = 1, . . . , k, i.e. only at the global minimum {U = WD{I © Xq^D"^ : 
Xq e U(no)}. So, there are no local traps in the kinematic landscape, and the remaining HiLil"-' + 1) — 2 
critical submanifolds are all saddles. 



3 Kinematic Critical Point Analysis of Phase-Invariant Landscapes 

Now consider the phase-invariant kinematic landscape Jp{U) | Tr(j4A^W^^^[/)p. Since Jp is phase- 

invariant, there exists a well-defined function Jp : PU(A^) K such that Jp ~ Jp o tt, where PU(A^) is the 
projective unitary group U(A^)/U(1)I and tt : U(A^) — ?► PU(A^) is the natural projection map depicted in 
Figure[T] One can then lift Jp to Jp : SU(A^) — M on SU(iV), the universal covering space of PU(iV), by the 
covering map p : SU(A^) — !• PU(7V). Here, Jp is just the restriction to SU(A^) of Jp, and it is easily verified 
that grad Jp([/) — gradJp(f7) and Hessj^ jy = Hessj^_[/ |T^^gu(Ar) for U e SU(A^). Many components of 
the critical set of Jp will admit exactly N disconnected images in the critical set of Jp. Choosing one 
representative such image for each connected component of Jp permits simplifying the analysis below. 
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U(l)^U(iV) 



SU(iV) ( ^"""'"''^^ z/NZ 





PU(7V) 



Figure 1: U(A^) and SU(A^) as fibre bundles over PU(iV). 



3.1 Distance Metric on PU(iV) 

There are various ways of deriving a phase- invariant landscape like Jp from one that is phase-dependent 
like Jf- a simple approach is to observe that min^gu(i) \\{U - (l)W)Af = 2\\A\\'^ - 2\Tv{AA'<W'^U)\. 
This provides a means to define a quotient metric on the projective unitary group PU(iV) from the metric 
d{U, W) = II (t/ - W)A\\'^ on U(A^). 

Another approach involves the adjoint representation of the unitary group, Ad : U(A^) — > Aut (u(A^)) C 
GL(u(Af)) ^ GL(Af2. jg gjygj^ Ad(£)A = UAU^ for any A € u{N), and where Aut (u(iV)) is the 
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2o| . With u{N) given the Hilbert-Schmidt inner product, 



group of Lie algebra automorphisms on u(A'^) 
{Ad{U)A,Ad{U)B) = {UAU\UBW) = (A, B), so that for each C/ G U(iV), Ad(C/) is an orthogonal operator 
on u{N). Furthermore, if W and U differ only by a global phase, i.e. W = e^^U, then Ad(U) ~ Ad(W^). 
Moreover, the kernel of Ad, i.e. Ad^^id), is the center of U(iV) 6, Cor. 5.2, pg. 129][2fl, Thm. 3.50] 
which is Z(U(iV)) = {e'^I} . Then Ad(C/) = Ad(W^) if and only if = e'^t/, so the image of Ad is 
a faithful representation of PU(A^) and Ad may be thought of as playing a role similar to the projection 
TT : U(A^) ^ PU(iV). 

Consider some target W E U(iV) and some B E GL(u(A^)) and define 



J({7)- II (Ad({/)-Ad(W-))oB||l,5, 



(11) 



where || • ||_ffs is the Hilbert-Schmidt norm on End(u(iV)) = , the space of all linear operators 

on u{N). Let {i^j} be the othonormal basis for u(A^) (with the Hilbert-Schmidt inner product), with the 
elements i\p){p\ for p = l,...,iV; ■^{\p){q\ + \q){p\) for 1 < p < q < iV; and ■j^{\p){q\ - \q){p\) for 
1 < p < q < N, where {\p)} is the canonical basis for . Then, 



J{U) = 2 Tr(B* o B) - 2 Tr (B* o Ad{W)* o Ad{U) o B) 
= 2Tt{B* o B) - 2^{B{nj),Ad{W^U)B{nj)). 



(12a) 
(12b) 



Let A E GL{C^) be an arbitrary A^xiV matrix, and let B be defined by B{n) = AQA^ . Then (B(rjj), Ad{W^U)B{nj)) 
can be computed for each of the basis elements as in Table [T] We find that 
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{B{nj),Ad(w^u)B(nj)) 


i\p){p\ 


\{p\A^W'<UA\p)\'^ 


^{\p){q\ + \<l){p\) 


mp\A''W''UA\p){q\A'<U^WA\q) + {p\A'fW^UA\q){p\A'<U''WA\q)) 


^,{\p){q\-\q}{p\) 


mp\AWUA\p){q\A''WWA\q) - {p\AWUA\q){p\A''WWA\q)) 


Table 1: Value of the diagonal element {B{flj),Ad{W^U)B{Clj)) for each type of basis element Clj G 


n{N). 





J{U) = 2 Tr{B*B) - 2 ^(B(^), Ad(f/t (13a) 

j 

N 

= 2Tr{B*B) -2^\{p\A^W^UA\p)f 

p=i 

-2Y^mp\A'<W''UA\p){q\A'<U''WA\q) + {p\A'<W''UA\q){p\A'<WWA\q)) 

p<q 

-2j2mP\^''W^UA\p}{q\A^U^WA\q) ~ {p\A^W^UA\q){p\A'<U''WA\q)) (13b) 

p<q 

N 

= 2TiiB*B) -2j2\{p\^^W^UA\p)\^ - 2j2{p\A^W^UA\p){q\A^U^WA\q) (13c) 

p=l pjiq 

N N 

= 2Tt{B*B) - 2^^(p|Atwt[/^|p)((^|^t[/t^^|(^) (13d) 

p=l q=l 

= 2^(B(0,),B(n,)) - 2| Tr(Atwt[;^)|2 (13e) 

j 

= 2\\A\\%s - 2| Tr{AA^W^U)\' = 2Jp(U). (13f) 

Consequently, the kinematic landscape Jp{U) = — | Tr(AA^ VK^J7)p on U(A^) is equivalent to the 

weighted Hilbert- Schmidt distance function on the subgroup of SO(A^^) given by Im(Ad) ~ PU(A''). 

3.2 Critical Point Identification 

Now, the differential of Jp at U G V{N), duJp ; Tf/U(iV) R is given by 

duJp{6U) = - Ti{AA^W^6U) Tr{U^WAA^) - Tr{AA'<W''U) Ti{6U^WAA'<) (14a) 

= ( Tr{U''WAA'')UAA'<W''U - Tr{AA^W^U)WAA'' ,6U) (14b) 

so that 

grad Jp(?7) = Tr{U^WAA^)UAA^W^U - Tr{AA^W^U)WAA^ . (15) 



8 



For Jp{U) < Pll^, Tt{WWAA'<) ^ 0, so grad Jp(J7) = if and only if AA'' Z = Z'<AA\ where Z = 
I T^^c/tM/AAt)! e U(iV). Using the notation from the previous section, D e U(iV) is such that AA^ D = 

with il'jj — ujj and ujf > ■ ■ ■ > ujj^; additionally, Cjf > ■ ■ ■ > oj"^ > Coq — denote the distinct eigenvalues 
of AA^ , with {rti, . . . ,nK,no} being the multiplicities of these eigenvalues. Then AA^ Z — Z^ AA'^ can be 
rewritten Vi^Z ^ Z'^Vi^ where Z = D^ZD. This same condition was considered in section [2] (and [8|), 
where it was shown to imply that [Z,AA^] and Z = f Af"f ffi Zq where t e G = U(ni) © • • • ® V{n^), 
A = Ai © • • • ® A^, A, = (-1,,,) ©I„._^, e U(n,), < i^^ < n„ and Zq G U„„. 

Drawing on the material above, we find that any critical point U of Jp with Jp{U) < can be 

written as U = | 'Ti'(uf\VAA^\ \ {TAT^ © Zo)D^ . This characterization is complicated by the presence of 
U on both sides of the equation, especially with regard to the phase factor on the right hand side. We may 
write U = e'''WD{fAr'< © Zo)D'f for any 6i e M, f e ^, A = -I^, © I„,_^, © • • • © -I^^ © I„,-^,, and 
Zq e U(no). Then [U^W,AA^ = [W'fU^AA^] = and WWAA'' = (AA^W^U)^ = e-^^^AA^W^U, so that 
gradJp(U) = 0. Hence, every such ?7 is a critical point of Jp, and they comprise connected critical sets 

Ca := {U = e^'^WDiTAV^ © Zo)D'' : e''^ € U(l),f e G, and Zq G U(71o)} (16) 
= {U = WD[{e'^rAr^) (B Zo]D^ : e**^ e U(l), f e and G U(no)} (17) 

for A = -I^^©I„i_^i©- • •©-It.^©I„^_y^ suchthat Tr(ri2(A©I„„)) > (hence Jp(t/) < WAW^). Such a critical 
submanifold Ca has the critical value Jp{U) = Tiin^ - Tr {Q^iA © I„J)^ = ( J^'^^^ n,Cj}f - ( E»"=i("» - 
2vi)Cjff = 4(X;r=i '^»w2)(X;r=i("i - for every [/ e Ca. Consider the set R := {e^^f Aft : e 

U(l) and f e ^} C U(iV-no) for this same fixed A. An element of i? has determinant e'(^""")''(-l)^ , so 
for it to lie in i?nSU(iV-no), we must have that 6 = 7r(2fc+X; ^i) /(N-no) for some k. Thus, RnS\J{N-no) 
is fractured into N — uq distinct connected components, each one being the image of the projection of R 
into PU(A^ - no). Let 9 be any of the N - uq admissible phases above. Then R := {e'^FAft : f e is 
one of the connected components of i?n SU(A^ — rto), and R — tt~^{p{R)), wherep and tt are the projections 
from Figure [TJ The projection p is a covering map and therefore is necessarily also a local diffeomorphism, 
and since it is a bijection on R, it is a diffeomorphism between R and p{R). The fact that i? is a dimension 
d = J2i=i {^1 ^ ^1 ^ i^i ^ ^iY) = 2 X]i=i ~ ^i) manifold diffeomorphic to 

U(ni)©...©U(nJ -Gr.,(C"0®---®Gr..(C"i (18) 



U(i/i) © U(ni - j/i) © • • • © U(z/«) © U(n« - v^) 



implies the same about p{R). Then, since the projection tt : \]{N — tiq) — > PU(Af — rtg) is a surjective 
submersion, R = 7r^^(p(i?)) is a dimension 1 + 2 X)i=i '^il"-* ~ '^O submanifold of U(A^ — uq). Moreover, R is 
a U(l) principal fibre bundle over Gr^.^ (C"i) © • • • © Gr^^ (C"") . Then R © U(no) is a dimension 1 + + 
2 X]r=i '^^("■i ~ ^i) submanifold of U(A^), and since conjugation and left translation are both diffeomorphisms 
on U(A^), Ca = R® U(no) is also a dimension 1 + + 2 X^JLi ~ ^0 submanifold of \]{N). 
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3.3 Hessian Analysis 

Given the form of the gradient of Jp in P3)) . by again extending the gradient vector field to C^^^ and 
differentiating, it is found that 

dc/grad Jp(JC/) = Tt{5U''WAA'')UAA''W^U + Tr{U^WAA'')6UAA''W''U 

+ Tt{U''WAA'')UAA'<W''SU ~ Tt{AA''W^U)WAA'<, (19) 

whence by projection onto the tangent bundle of U(A^), 

Ressj^,u{SU) ^ Vsu grad Jp = - Ty{AA^ WUU)W AA^ + Tr{5U^W AA^)U AA^W^U 

+ i| Tv[U^WAA^)5UAA^W'^U + Tv[AA^W^U)WAA^ SUU'^ 

+ TriU'^W AA'^)U AA'^WHU + Tr(AA^T4^^[/)t/U[/VKAA^|. (20) 

On one of the critical submanifolds Ca, the Hessian is given by 

Hessj^,(7((5J7) = Tt:{AA^W^U)5UU^WAA^ +Tt:{AA^W^U)WAA^uUu 

^21:y{AA'^W^U)WAA\ (21) 

Since U = e^'^WDiVKV'^ ® Zo)D^ can also be written U = {-e''^)W D{T{-k)t^ ® -Zq)D\ we can choose 
to always write such a critical U such that Tr(n^(A Ino)) > 0- We seek an eigenvalue decomposition of 
Hessjj,_(7, and find that '^SU = Hessjp_c/((5C/) if and only if 

7^ = Ti{n^k){Yifk + ll^Ar) - 2 Tr(172Ar)17^A (22) 

where y = (ft ® I„ jDt^/t^c/ £)(f ® I„o) e u(7V) and A = A I„„. We begin by looking at the nuh 
space of Hessc/„/p, where 2TT:{n'^ kY)VL^ K = Tx{n^k)[YVL^k + ^^kv). Then Tr(rj2A)(A,wf + Ajw|)ry = 
2Tr(f2^AF)AiWj^%. For i ^ j, and since Tr(ri^A) 7^ by assumption (else Jp{U) = \\A\\^), this implies 
that either X^ujf + Xjui'^ = or Y^j = 0. For i = j, we get TT{n'^k)ujfYii = Tr(f22Ar)a;f , implying that the 
diagonal elements of Y are all identical where ujf 0. So Y = (ij/dljv-no) ® 0no + Y, where j/d G K and F 
is of the block diagonal form shown in ([9]). Hence the nullity of Hess ^(7 is Ao — 1 + + 2^Vi{ni — vi). 

For the remaining eigenvalues 7 7^ 0, we have 7^ - i:Y{V?k){YVL^k + O^AF) = -2TY{n'^ kY)n'^ k 
or ^Yij — k){\iLjjf + XjU}^)Yij = —2TT{fl'^kY)ujfXiSij. Since the right hand side is diagonal, when 

i ^ j either 7 = Tr(ri2A)(Aiwf + AjCj^)^ or Yij = 0. Since we chose A such that Tiift'^k) > 0, 7 < if 
XiUjf + XjUjj < for i ^ j, i.e. whenever 

• Xi = Xj = —1, which happens for (X^jLi ^i) (SJLi ~ ^) ordered pairs (i, j), 

• Xi — —1, Aj = +1 and ivf > ujj, which happens for X]fc>i>i ^k{ni ~ vi) + rio X]r=i ordered pairs (i, j), 
or 

• Xi= +1, Aj = -1 and < uj'j, which happens for I]i<fc<i("fc " '^fc)t^i + "0 Ya=i = J2k>i>i ^kim - 
vi) + TiQ X^jLi ordered pairs (i, j). 
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Likewise, 7 > if AjW? + XjOJj < for i ^ j, i.e. whenever 

• Aj = Aj = +1 and at most one of ojf and 00 j is zero, which happens for 
{N-no- Er=i ^0 - no - 1 - X^f^i vi) + 2no (A'' - no - Ylk=i ^k) ordered pairs 

• A, = —1, Xj = +1 and cof < ojj, which happens for Yli<k<i — ui) ordered pairs or 

• A, = +1, Xj = -1 and cjf > w?, which happens for Ylk>i>i('"'k - Vk)n = Yli<k<i ^kini - n) ordered 
pairs (i, j). 

On the other hand, when i = j, we get 7Fii-2 Tv{n^A)XiUjfYii = -2Tr{Q'^AY)XiUjl If 7 = 2 TT{n'^A)XjUJ^ = 
2Ti{n'^ A)XkUjl, andy = ^ - \k){k\), then Tr(f]2AF) = and -fYjj +2 Tr{n^A)XjOj]Yjj = 2Tr(n^AY)Xjw]. 

This accounts for an additional X^iLi n:ifix(t'i — 1,0) = ^^-i Vi — #{i : 1 < J < k and > 0} eigen- 
vectors for 7 < and Yli=i niax(ni — z/^ — 1,0) = N — tiq — Yli=i ~ #{i • 1 < j < and Vj < 
rij} eigenvectors for 7 > 0. Finally, consider the case where Y is diagonal and Tr(f2^AF) ^ 0. Then 
(7l-2Tr(n2A)02A)F = -2 Tr(02 Ay)n2^, so for 7 ^ -2 Tr(02 A)Aja;2 for all j, Y = ^^^gg'/g^^^ . This 

yields Tr(n2Ay) = Tr(f22Ay) Tr ( 2Tr(o4"o^A-7i[ ) ^"'^ ^^"^^ ^® assume that Tr(02Ay) ^ 0, it follows that 
7 must satisfy 



/(7) := Tr . = > = 1- (23) 

\^2Tr(r!2A)l]2A-7iy ^ 2 Tr(f22A)A,w2 _ ^ 

If there are I distinct values ryi < • • • < 77; of 2 Tr(f22A)Aja;2 with w? ^ 0, the function / has / solutions 
7(7) = 1: 71 < ryi and 77^-1 < jk < r/k ioi k = 2, . . . ,1. Note that 7/; = is one of these solutions, but should 

be neglected since we have already characterized the null space of Hess jp^u- Letting Y = 2'af^A 

where a > is a normalizing factor, we see that Tr(02Ay) = ia ^ 0, Tr(y) = 0, and that y is a 
diagonal eigenvector of Hessjj,_(7 with eigenvalue jk- The number of eigenvectors with negative eigenvalue 
contributed in this way is equal to #{j : 1 < j < k and i^i > 0}, while the number of eigenvectors with 
positive eigenvalue is equal to #{j : 1 < j < k and fj < rij} — 1. 

To summarize, away from the global maximum, the function Jp has connected critical submanifolds 

C'a := {U = e*''T4^£'(fAf^ ® Zo)D'' : e''^ e U(l),f e Q, and Zq G U(no)} (24) 

for A = ® Ini-1.1 © • • • © -Ii., © such that Tr{n'^A ® I„J > (hence Jp{U) < \\A\\^). Such a 

critical submanifold has dimension c^a = 1 + + 2 ^ i^^^-^ {rii — Vi), which is identical to the nullity A/q of the 
Hessian on Ca- Tallying up the dimensions of the negative and positive Hessian eigenspaces, we find A/1 = 
(Er=i ^0^ + 2 Efc>(>i {m + 2no X;r=i '^i and Af+ = [N - no - uif + 2no {N - no - X^Li ^fc) + 
2Si<fc<i '^k{ni — vi) — 1. As a result, M- = if and only li Vi = Q for all i = 1, . . . , k, i.e. only at the 
global minimum {U = e*^VF£'(I© Zo)D^ : Zo € U(no) and G [0, 27r)}. Furthermore, Af+ can be zero 
for such a critical point only if uif > X]i=2'^?' "o = 0, and A = 1 © — I^r-i. The critical value of such a 
point is \\A\\^ — {ujf — Yli=2 ^?)^ < ll^ll^' in this very special circumstance, there exists a local maximum 
submanifold distinct from the global maximum set {U : Jp{U) = ||j4||^}. Since the intended use of this 
function is for minimization, such a local maximum submanifold does not act as a "trap" . 



11 



3.4 Global Maximum Set 

Finally, consider the global maximum set {U : Jp{U) = \\A\\^} = {U : Tr(AAtH^t[/) = o}. This set does 
not admit analysis by the methods used thus far, so a different approach is required. Let F : 1LJ(A'') M? be 
given by Fi := ^Tv{AA^WW) and F2 := 3 Tr(AAtVFt{7). Then, 

duFi{SU) = ^Tr{AA^W^SU) = ^{WAA^ - UAA^W^U,SU) (25a) 

duF2{SU) =QTi{AA^W^6U) = ^{iWAA^ + iUAA^W^U,dU) (25b) 
so that the gradients are given by 

gTadFi{U) = ^{WAA^ -UAA^W^U) (26a) 

gradF2(t/) = ^{WAA^ + UAA^W'^U). (26b) 

Thus, dF is surjective except where there exists (a,/?) ^0 gM? such that Q!gradFi(f/) = ^gradF2(J7), i.e. 
where 

e^t'U^WAA^ = e-"^AA^W^U, (27) 

where e"^ — {a—i/3)/\a+if3\. As wo have already seen, this equation implies that U ^ e"t>WD{rAP®Zo)D^. 
For such a U, Ti{AA^W^U) = e'"^ Ti{n'^A), so the only way for Ti{AA^W^U) to be zero is if n'^ is orthogonal 
to one of the possible A's. This means that for a given A'' and no, would have to lie in the union of the 
2JV-no-i hyperplanes which form the orthogonal spaces of the A's within the space of real diagonal matrices. 

For a given and no, the collection of all A's for which dF is surjective at F{U) = [hence {U : F{U) = 0} 
forms a codimcnsion 2 submanifold of U(A^)] corresponds to an open dense subset of C'^^-^ . Now, at a point 
U such that Tt{AAHV^U) = 0, 

Hessjp,cj(<5[/) = - Tr{AA^W^SU)WAA^ + Tr{SU^WAA^)UAA^W''U (28a) 
= -[^Tt{AA^W^SU) +mTT{-iAA^W^SU)]WAA'< 

+ [^Ti{6U''WAA'<)+mTi{-idU''WAA'<)]UAA^W^U (28b) 

= -^[(WA^t - UAA^W^U,SU) + i(iWAA^ + iUAA^W^U,SU)]WAA^ 

+ ^[{WAA^ - UAA^W^U,5U) - i{iWAA^ + iUAA^W^U,SU)]UAA^W^U (28c) 

= -i|(iy^A^ - UAA^W^U,SU)(WAA^ - UAA^W^U) 

+ {iWAA'' + iUAA'^ W'^ U, 5U) {iWAA"^ + iUAA'' W"^ U) } (28d) 

so that the Hessian is rank 2 except where there exists (a, /3) 7^ e such that q:(W^^^ — UAA'^W^U) = 
[}{iWAA^ + iUAA'^W'^U). This is exactly the condition just considered for the surjectivity of dF, so that 
the Hessian is rank 2 if and only if dF is surjective. So, on the open dense set of A's for which this happens. 
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the maximum set is a nondegenerate (in the Morse-Bott sense), codimension 2 submanifold of U(A^). Since 
the other critical points also comprise nondegenerate submanifolds, we conclude that for these A's, Jp is 
a Morse-Bott function. In the particular case that AA^ is fully degenerate (e.g., A = 1), it is found that 
the maximum set of Jp is a nondegenerate submanifold if and only if N is odd. However, when N is even, 
arbitarily small perturbations of A about I are sufficient to obtain a Morse-Bott function. 



4 Dynamical Critical Point Analysis 

Now that we have elucidated the structure of the critical sets of the kinematic landscapes Jp and Jp, 
we return to the problem of characterizing the critical set of the dynamical landscapes J = J o Vp- Let 
M C U(iV) be one of the critical submanifolds identified in the previous sections. It can be proved (see 
Appendix El that Vr : K -)■ U(A^) is C°°. In addition, since U(A^) is finite-dimensional, if Vp{E) £ M 
then (dfVr)"^ (Ty^(£)A^) has finite codimension, so is closed and has a closed complement (i.e., it "splits"). 
Therefore, away from singular points of Vp (i.e., those £" € K such that dgVp is rank-deficient), Vp is 
transversal to Jv[ and by the transversal mapping theorem Vrj^^{M) is a Hilbert submanifold of IK, 
TeiVp^^M)) = {deVT)-^{Tv^^£)M), and codim (l/f ^(X)) = codim(X). 

Let £ S K be a regular critical point of J, i.e. such that grad J{£) = and dgVp is full rank. It may be 
seen that at such a point, the Hessian of J is given by Hessj^ = (dsVp)* o Hessj_yj,(f:) o(d£VT)- Let As be 
the linear operator on Ty^(^£)\J{N) given by 

As = {dsVp o {dsVpyy o Hessj,y^(£) o{dsVp o {dsVpYy . (29) 

Since dsVp is assumed to have full rank, we may apply Sylvester's law of inertia [o^ to conclude that 

As and Hessj v'T(£) have the same numbers of positive, negative, and zero eigenvalues. Let {ir]j,Vj)} for 

_ i_ 

j = 1, . . . ,iV^ be the eigenvalues and eigenvectors of As, and let Zj — (dsVp)* o (^dsVp o (dsVp)*) ^Vj. 
Then 

Hessj^ Z, = {dsVp)* o Hessj,y^(£) o{dsVT o (dsVpYy^j (30a) 
= VjidsVp)* o {dsVp o (dsVpyy^V, = T]jZj (30b) 

so that {{rij, Zj)} for j — 1, . . . , N'^ are eigenvalues and eigenvectors of Hess j ^. Since Hess j ^ is self-adjoint, 
any other eigenvector Z must be orthogonal to the {Zj}. Also, note that since the {V,} span Ty^(^g)\]{N), 
the {Zj} span Range {{dsVp)*). Then, for any X e Ty^(£)U(iV), = (Z, {dsVp)*{X)) = {dsVp{Z),X), so 
that dsVp{Z) = and therefore Hessj^ Z = 0. Thus, Hessj^ has infinitely many eigenvalues; N"^ of them 
are identical to the eigenvalues of As, and the remaining infinite number of eigenvalues are all zero. Since 
J has no local traps, we can conclude that J has no local traps among the regular critical points. From 
the transversal mapping theorem we find that T£(V^^(A^)) = (d^Vr)"^ (Tyy(£)A^) , implying that for any 
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/ G K, wc have 

Hessj£(/) = ideVT)* o Hessj,vv(£) "(deVrKf) = ^ dsYrif) G kerHessj,^ (31a) 

^d£VTi.f)eTv^^£)M (31b) 

^ / e (d£yT)-^(Tv,(£)A^) (31c) 

^/eT£(vv'(X)). (3id) 

Hence, the null space of Hessj^ is identical to Tf (V^j7^(A^)) , the tangent space to the critical submanifold. 

In the case where the Hamiltonian takes the dipole form H{t) = Hq — £{t)iJ. for any £ E L^(IR+;R), 
the Frechet derivative of Vt is given by dsVriSS) = tVV(f) Vt\£)nVt{£)S£{t)dt 5]. Then the adjoint 
operator of the derivative is d£V^{A){t) = -^Tt {A'fVT{£)Vt\£)^J.Vt{£)) for any A £ Ty^(£)U(iV), and 
the operator norm of this adjoint is uniformly bounded by jjdflj^ll < -s/Tll/iH- For any smooth 'kinematic' 
function g : U(A^) — > M, let g = 5 o Vr be the corresponding 'dynamical' function on L^(R+;R). Then 
gradg(f) = df^^^ (grad <7(l^T(f ))) and |1 grad5(£:)|| < \\d£Vf\\\\ gve.dg{VT{£m < VT||Ai|| || gradg(yT(f ))|| • 
Since g is smooth, || grad^jj is continuous over U(A^), so that since U(A^) is compact, || gradgH is uniformly 
bounded. Therefore, || grad^jj is uniformly bounded over i^(R+;R). For any dynamical quantum control 
landscape constructed in this way, the slope of the landscape is uniformly bounded by some constant. 



5 Landscapes Based on Intrinsic (Geodesic) Distance 

The kinematic landscapes Jp and Jp considered above are based on the Euclidean (or norm) distance on 
U(iV) and PU(A^), respectively. We now describe two additional distance measures based on the intrinsic 
distance between operators in U(7V) and PU(iV) under the Riemannian metric induced by the real Hilbert- 
Schmidt inner product on C^^^. 

The first of these distance measures is quite simple to define. Since the chosen Riemannian metric is 
bi-invariant on U(7V), any geodesic starting at C/ G U(A^) is of the form j{s) — Ue^'^ for some A G u(A^). 
To find a geodesic joining U to some target W G U(iV), let C/e^ = 7(1) = W, so that e^^ = WW and 
A — \og{WW). This matrix logarithm is not uniquely defined, but the length of the geodesic 7 defined on 
the interval [0, 1] is given by £[7] = Il7'('5)|| ds = ||^||. The minimum such length is obtained by taking 
A = \og{WW) from the principal branch of the logarithm so that all eigenvalues lie in (— itt, itt]. We then 
define the landscape as Jg{U) :— ^|| log(C/^VK)|p. Then the gradient of Jq is given by (see Appendix [Bjl 

grad Jg(?7) = -UlogiU^^W). (32) 

As most numerical matrix logarithm routines (e.g., the logm function in MATLAB) compute the principal 
branch, they provide a ready means to obtain both the landscape value and the gradient. Since the norm 
of grad Jg is the distance to the target, this vector field is only zero at the target, i.e. the global minimum 
of the landscape. Hence, there are no traps or saddles. The gradient field has the property that it is 
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discontinuous and multiply defined at the cut loci of U(7V) (where the spectrum of U^W contains —1), but 
this is not a problem for an optimal control algorithm since the matrix logarithm routine will have to choose 
one from among the possible solutions, all of which describe minimal geodesies to the target that are equally 
satisfactory. 

A phase-invariant version of Jg may be constructed analogously by considering minimal geodesies on the 
projective unitary group PU(iV) ~ U(A^)/U(1), or equivalently by defining Jgp{U) :— min^ i|| \og{e^'^U^W)\\'^ 
on U(-/V). It may be shown that 



Jgp 



1 



mm — 

keiN 2 



log dei{U^W)~Tru'^W^ 
grad Jgp(C/) ^ -C/jlog {e^ dei{U^W)-TrU^W^ 



N 



Tr 



los 



(33a) 



(33b) 



where k in (|33bp is the minimizer from ([33a|). With this minimizing fc, the trace in (j33b|) will be zero, so 
that 

grad Jgp(?7) = -U\og [e^ det(L/^l^)"^J7^VF) . (34) 

As with Jg, the norm of grad Jgp is the distance to the target, and this vector field is only zero at the 
target, i.e. the global minimum of the landscape. Hence, there are no traps or saddles. One downside to 
this landscape is that it appears that all N possible values of k must be tried in order to find the minimizer 
of p3al) . This behavior has a topological interpretation on PU(A^). Since the fundamental group of PU(iV) 
is 7ri(PU(7V)) = Zjv = Z/A^Z, there are exactly A^ homotopy classes of paths connecting 7r([/) to the target 
n{W). Within each of these classes is a unique minimal geodesic, and these A^ minimal geodesies are identified 
by the vectors C/|log(e^det([/tVl/)-^J7^H/) -^Tr \og{e^ deiiU'^WyTrWW^ l| indexed by A: . 

A distance metric based on intrinsic distance could in principle be applied to the case where only some of 
the states are important, analogous to Jp and Jp where A is rank deficient (e.g. where A is a projector). This 
is equivalent to computing the geodesic distance between points on the Stiefel manifold Vfe(C''^) ~ \]{N)/lk(B 
lJ{N — k) or on its projective cousin T4(C^)/U(1). However, the two-point geodesies on these spaces are non- 
trivial to compute. The problem requires solution of a boundary value problem or an optimization problem 
to find each minimal geodesic. For that reason, these intrinsic distance metrics may not be practical for this 



scenario. 



6 Summary 

This work presented an expanded analysis of landscapes Jp and Jp, which are based on the Euclidean 
distances between unitary operators in U(A^) and PU(A'^), respectively. The expansion appears in several 
ways. First, additional freedom has been allowed in the landscape functions themselves, by admitting A 
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matrices that are rank-deficient. Landscapes based on these rank-deficient A matrices measure the distance 
between unitary operators by their action on a subspace of the full state space. This can be the desired 
objective for designing a quantum information processor, for example, where only this subspace of the state 
space is to be used for the quantum register. This additional freedom in defining the landscape is consistent 
with the principal finding of earlier work on landscapes of this form: they have no suboptimal minima 
(i.e., "traps") that could impede a deterministic optimal control algorithm (such as gradient descent) from 
reaching the global minimum. 

In addition to broadening the families of landscapes for consideration, we have provided more detail on 
the structure of the critical sets and the behavior of the landscape functions at these critical sets. The critical 
sets were shown to generally be disjoint unions of critical submanifolds and we have described the structure 
of these submanifolds, either as products of Grassmannians and a unitary group in the case of Jp, or as a 
U(l) principal fibre bundle over such a product in the case of Jp. Furthermore, we have shown that these 
critical submanifolds are generally nondegenerate in the Morse-Bott sense, so that the kinematic landscapes 
are generally Morse-Bott functions. 

These results were related back to the corresponding dynamical landscapes through the control-to- 
propagator map Vt, implicitly defined by the Schrodinger equation, that takes a control function as input 
and returns the final time unitary evolution operator. This map was shown to be infinitely Frechet dif- 
ferentiable, leading to the conclusion that, away from the singular points of Vt, the level sets and critical 
sets of the dynamical landscapes are C°° smooth, finite codimension submanfolds of the infinite dimensional 
control space K = L^(IR+;IR). Also, the number of positive and negative Hessian eigenvalues (and therefore 
the characterization as a minimum, maximum, or saddle) was shown to be identical for a kinematic critical 
point and a regular point of Vt that maps to it. This behavior implies that no traps exist in the dynamical 
landscape among the set of regular points of Vt- Furthermore, Morse-Bott nondegeneracy of the critical set 
is also preserved away from singular points of Vt, which can be important for certain numerical landscape 
exploration methods such as second order D-MORPH 

Finally, two additional landscapes were introduced that are based on the intrinsic or geodesic distance 
between operators in IJ{N) and PU(iV), respectively, rather than Euclidean distance. These kinematic 
landscapes have the desirable property of having no critical points except at the target. These landscapes 
may allow for more efficient performance of optimal control algorithms over Jp and Jp, since the latter have 
many saddle points where the gradient is zero. 
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A Differentiability of U{T, 0) With Respect to the Control 



Let M(A^) C C^^^ denote the space oi N x N Hermitian matrices endowed with the real Hilbert-Schmidt 
inner product {A,B)hs = ^TriA'fB), and let H(7V) = L'^{R+;M{N)) denote the space of all iV x iV 
square-integrable time-dependent Hamiltonians with inner product {Hi,H2)l-- = Jf^^TiiHl{t)H2{t)]dt. 
Let Zt '■ M{N) — >■ U(A'') be the map, defined implicitly through the Schrodinger equation, that takes a time- 
dependent Hamiltonian H(-) G 'E.(N) and produces the corresponding unitary time-evolution operator at 
time T: U{T, 0) e V{N). This map is well-defined over the entire domain because of the absolute convergence 
of the Dyson series over H(iV): 

Zt{H)=I+(-1] f dhH{ti)+ (~^] f AhH{h) r dt2H{t2) + .... (35) 







In this appendix, we will prove that Zt is well-defined and infinitely Frechet differentiable over H(Ar). A 
corollary is that the map Vr : K U(iV) defined in the body of the paper is infinitely Frechet differentiable 
over all of K = L2(k^. k). 

Lemma 1. /// : [a,b] M. is integrable, then 

t t / b \ ^ 

dti fih) f ' dt2 /(<2) •••/'"' dt„ /(i„) = ^ I / .fit) I . (36) 



Proof. Note first that (l36l) holds trivially for n = 1. Suppose that it holds for n = m. Then 



dhfih) / dt2f{t2)' 



dtm+1 f{tm+l) 



1 

m 
1 



— / fih) ( / fit)dt] dh 



m\ dti 



1 



1 



m + 1 
b 



f{t)dt 

m+1 



dti 



(m + 1)! yj, 

and the lemma follows for arbitrary n G N by induction. 



f{t)dt 



(37a) 
(37b) 

(37c) 
□ 



Definition 1. For integrable matrix functions Ai : [0,T] — > C^^^ and for integrable real-valued functions 
at : [0,r] — >■ M, we will use the following short-hand notation for the Dyson-esque terms 



Tt[^i,A2,...,A„] := / diiAi(ti) / dt2^2(i: 



dtnAnitn) and 



TT[ai,a2, . . . ,a„] := / diiai(ti) / dt2a2{t2) 



dtn 0>n (^n) ■ 



(38a) 
(38b) 
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Lemma 2. If gi,- ■ ■ ,gm € Z/^(IR+;K+) are non-negative square-integrable functions on [0,oo), then 

'TT[gi,g2,---gm]<T"'^^gi\\L^---\\gm\\L^ (39) 

Proof. Since the gi's are non-negative functions, wc get the inequahties 

'^T[gi,g2,---gm]= dtigi{ti) dt2g2{t2)--- dtmgm{tm) (40a) 
Jo Jo Jo 

< [ dtigiiti) [ dt2 52(i2)--- / dt„,g^{t^) (40b) 
Jo Jo Jo 

<T"^/^(^£ dt.giit^)^ ^£dt2gl{t2)^ dtmglitm)^ (40c) 

<T"^^'\\gi\\L^---\\gm\\L^ (40d) 

by extension of the integrals out to the interval [0,T], followed by application of the Cauchy-Schwarz in- 
equality, and finally extension out to [0, oo). □ 

Lemma 3. If f,gi, . . . ,gm S -L^(]R+;IR+) are non-negative square-integrable functions on [0, oo), then 

'^T Ij^-jjjj 91 , f^-jjj 92,f,---J,gm, f^-jj) - Yl"" ^i ll-^ll^" "^T[gi,g2,---gm] (41) 

terms terms 13m terms 

Proof. Let ai = (5o + ■■■ + (3i + i -\- 1 for i = 0, ... ,m. Then, using Fubini's theorem, we may rearrange the 
order of integration as follows: 

'^t [/, ••^■,/ , gi, f,..j,fj , 52,/, • • • , f,gm,L_-jJ) 

Po terms /3i terms Prn terms 

rT rt<To-2 /■*<^o-i /■*''o 

= / dtif{ti)-- - / dt„„-if{t„„-i) / dt„„gi{t„J / dt„„+if{taa+i)--- 
Jo Jo Jo Jo 

••• / dtai-i f{ta^-i) / dtaig2{tai) / dt^^+l f {ta^ + l) ■ ■ ■ 

Jo Jo Jo 

••• / dtam-1-1 f{ta,-i-l) / dtam-i 9m{tam-i)x 

Jo Jo 

X / '"''dt,„-i+i/(ia„-i+i)--- / " dt,„_i/(t,^_i) (42a) 
Jo Jo 

1 f'^ f^'O /■*<'m-2 

= rr"' I / dtao5i(i-7o) / dt„^ g2{t^,) ■ ■ ■ I di<,^_i 5„ (<<,„_ J x 
lli=0 "'0 Jo Jo 

X l'^^ /(t) dt^ y^''" fit) dt^ ■■■ fit) dt^ ^ (42b) 
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where the last step follows from Lemma [T] Then, since / is a non- negative function, we get the inequality 



^0 terms /3i terms (5m terms 

^ rrm^ ^ I f / /Wdt) / dhgiiti) [ dia .92 (i2) • • • / dt™5„(i„) (43a) 



Ti^'^' II iiE/3- 

< prm ^ , ||/||j,2 'TT[gi,ff2,-.-,gm] (43b) 
lli=0 Pi- 

by first extending the / integrals to the interval [0, T]. and then invoking the Cauchy-Schwarz inequality. □ 
Lemma 4. 

Proof. First, observe that for any a > 0, 



Then it follows that 



(45b) 



E n^ = EE- E E — ^r^:^. (*> 



EE-- E . . ' .; ^.-2 X. (46b) 



ci=oc2=o c,_2=o ci!c2! • • • c,— 2! (^p - ! 



P _ l)P-Ci 

E „ M (46c) 



(46d) 
□ 
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Definition 2. For m ^ 0,1,2, ... , let B'^(TM.{N); C^^^) denote the space of bounded m-multilinear oper- 
ators from (TH(iV))'" = TH(Ar) x TUiN) x • • • x TH(Af) to C^^^, with the norm 

{\\SHj\\jto} \\o^i\\---\\on„,\\ 
for each A e B™(TH(Ar); C^x^) . Then let $T,m : H ^ B™(TH(iV); C^x^) 6e defined by 



^T,m{H){5H^,5H^,--- ,5Hm):=Y, 

n=jn aoH \-am='n—^ n^Sm 



X 



X '^t[H,. . . ,H, 5H^ii\,H, . . . ,H, 6H„t2),H, . . . ,H,6H„/jn),H, . . . ,H] (48) 

ao terms ai terms am terms 

where Sm denotes the symmetric group on m elements (i.e., the group of permutations of m elements). For 
m,q = 0,1,2,..., let '^T,m,q : H © TH -i^ S"*(TH(iV); C^^JV-j defined by 



OO / 

^T,m.g{H,6H){6H,,6H2,--- ,6H^):= Y 12 E E( 

n=m-\-q ao + ---+a^_|_g \-bm.=Q TreSrr ^ 

n — m — g 



h 



I 



xrT[H,...,H,A,,H,...,H,A2,H,...,H,Am+q,H,...,H], (49) 

ao terms oi terms o,m+g terms 

where 

bo terms bi terms 

...,5H,5H^(^),5H,...,6H] (50) 

6m terms 

Lemma 5. $T,m awd ^T,m,g aJ^e well-defined since their defining sums converge absolutely, and for each 
if e H and 6H e TH, ^T,m{H) and ^T,m,q{H){SH) are bounded m-multilinear operators. 

Proof Let f{t) = \\H{t)\\, gi{t) = \\Ai{t)\\, hj{t) = \\5Hj{t)\\, and h{t) = \\5H{t)\\. Then 

WTtW, . . . ,H ,A^,H,. . . ,H,A2,H, . . . ,H,Am+q, H,. . . ,H]\\ 
^ V ' ^ V ' ^ V ' 

ao terms ai terms am+g terms 

<'^T[L_^^,gi,L_-2£,92,f,---,f,gm+q, (51a) 

ao terms ai terms o-m+g terms 

,52, • ■ • ,3r»+q] (51b) 

ni=0 

= ^^^7+5-^11^^11^2™ ^Tt [ fe, ■ ^- , /i , , /t, - ■ , ft ;, K(2) ,h,...,h, h^^^) , h,. . ^.,h ] (51c) 

^ ^ bo terms bi terms bm terms 
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< 7rT"^ ,, l|g||~^l|/^Hi.||feilUH|fe2|U^ • • • (51e) 

< ^„,+, 7rT"^ ^ ''||^HH|,H,5gi|U.||'^g2HL^---pg,n|U^ (51f) 

by appealing to Lemmas [5] and |31 Therefore, 
\\-^T,mAH, SH){5HuSH2,...,5H„,)\\ 

< y y y y — x 

X ||i7|r2^"-«||<5i/||^,||(5iJi|U2||5i72|U2 • • • ||<5i7„,|U. (52a) 
= E '"'^"^'";,;!'(' ! 2 -I)!a!''' ^ ||g||-'"-^H>5g||i.||<5gi|U.||'^g2|U^ ■ ■ ■ m,^W (52b) 

n—m-\-q ^ H)-H- 

m\T'^{m + lY f{m + q+l)VT 



■exp[ "'^ ' \ ^'- \\H\\l2 \\SH\\1,\\SH,U2\\SH2\\l- ■ ■ - WSH^U. (52c) 



so that the sum converges absolutely, '^T.m,q{H,5H) is a bounded m-multilinear operator for each H £W 
and 5H € TH, and 

\\^'T,mAH,SH)\\ < ^_A__^exp r- \ WHh^j WSHWl, < oo. (53) 

Since ^T,m{H) = 'i'T.m,o{H,6H), this conclusion also holds for ^T,m{H). □ 
Theorem 1. Zt is infinitely Frechet dijferentiable, i.e. C°° , everywhere on H. 

Proof. We begin by establishing that $T,m is Frechet differentiable for each ?7i = 0, 1, 2, . . . , with derivative 
*&T,m+i- Observe that 

$T,™+i iH){6Hu...,SHra,SH) = *T.™, 1 (H, 5H) (dH^ , . . . , 5if,„) (54) 

and that since the defining sum for ^T,m{H + 5H) converges absolutely, it may be rearranged as 

oo 

<^T,m{H + 5H){5Hi, 6Hm) = E '^T,m..'i{H, 5H){5Hi, ...,6H^). (55) 

9=0 
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Then by appealing to the bound of ^'r.m.q in (I52cp . we get 

||$T,,n(-ff + SH) - <i>T,m{H) - <^T,-m+l{H){- , ■,dH)\\ 

- sup \\<^>T.rr.{H + SH){SHi,...,dHr„)^'i>T.rn{H){SHi,...,dH„,) 

{\\SH,\\ = 1} 

- <i>T,m+liH)iSHu . ■ . ,SHmJH)\\ 

oo 

^ *T„„,g(if, SH){SHi,. . . , SHm) 



Hence, 



sup 

{\\SH,\\ = 1} 



q=2 

OO 



< sup y^\\^T,mAH,SH)iSH,,...,SH„ 

{II 5//, 11=1} 



< f; ""^ti;!^"' »p + „ 1 urn' 

q=2 



n<n+qq\ 



< 



E 



m!TT^ (m + 1)« ^ _^ / (m + q + 1) 
^m+g(g_ 2)! 



g=2 



■ exp 



\H\\j, \\5H\\% 



L2 



exp 



X exp 



{m+l)VT 



\H\\r2 X 



\SH\\l2 cxp\^\\H\\l2 



L2- 



lim 

11.5^11^0 



\m\\ 



= 



(56a) 
(56b) 

(56c) 

(56d) 

(56e) 



(56f) 



(57) 



and therefore $T,m is Frechet differentiable with derivative $T,m+i- Since Zt = ^t,o, this impUes that Zt 
is infinitely Frechet differentiable, and that the m'th derivative of Zt is ^T,m- □ 

Lemma 6. Let "H : IK — ;> ]HI(A^) be defined by 'H{£){t) — Hq — fJ-^it) for some fixed N x N Hermitian matrices 
Hq and /i. Then % is infinitely Frechet differentiable, i.e. C°° . 

Proof. Let C : K ^ 6(TK;TH(iV)) be defined by C,{£){5£){t) ^ -^fi6£{t). For each £ e K, C{£) is linear, 
and ||C(f)(<5f )|!l2 = llA^I|ffs||<5f IIl^, so ({£) is bounded. Now, 



\m£ + S£)-n{£)-a£){S£)\\ _^ 



\S£\ 



(58) 



so that C is the Frechet derivative of H. Since C, is constant (i.e., ({£) is the same linear operator regardless 
of which £ € K is input), the higher Frechet derivatives also exist and are all equal to zero. □ 



Theorem 2. Vr = Zt o "H : K — > V{N) is a composition of C°° maps and therefore is itself a C°° map. 
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B Gradient of Jq 

Several steps in Section [5] require differentiation of expressions involving the matrix logarithm. Since the 
expressions to be differentiated are all similar, this appendix will demonstrate the computation of the gradient 
of the kinematic landscape Jq, as the other variations follow along the same lines. To this end, we fix some 
target W G U(A^) and U E U(iV) and define the following maps: 

/ : u(iV) ^ M fiA) = iTr(AtA) = A) = 

g : U(iV) ^ u(A^) g{V) = log{V) (59) 

h : \J{N) U(7V) h{U) = WW. 

Then Jq — f o g o h. To compute grad Jg, we will determine the various pieces of the differential 

duJaiSU) = dg(h(u))f o d^^u^g o duh{6U) (60) 

and then recombine them as 

grad Jg(;7) -d^ro 4(^)3* ograd/(5(Mf/)))- (61) 

First, let us examine /(a) = ^{A,A). Because of the bilinearity and symmetry of the inner product, 
we can easily establish that the differential dAf '■ u(A^) — M is given by dAf{6A) = {A, 6 A) Hence, by the 
definition of the gradient, 

grad/(A)-A (62) 

For h{U) = U^W, it is easy to see that 

duh{SU) = dU^'W ^ -U''6UU''W. (63) 

Then to compute the adjoint of duh, we write for an arbitary SU G TijV{N) and 5V E Th(^if)lJ{N) 

{duh* {6V) , 6U) = {SV, duh{SU)) = {SV, -U^SUU^W) = { - USVW^U, SU). (64) 

Since 5V E T[;tH'U(iV), it can be written as SV = AWfW for some A E u{N), so that -USVW^U = -UA E 
T[/U(7V). Therefore duh* : Tft([/)U(iV) -> TuV{N) is given by 

duh*{SV) = -USVW^U. (65) 

Consider now the differentiation of g{V) = log(y), which needs to take into account that the N x N 
matrices form a non-commutative algebra. We begin by rewriting the definition of g as exp og — id, and 
differentiating both sides. Using the discussion of Kronecker sums in [3], we get 

SV^f exp{g(V)s) dvg{dV) exp{g{V){l- s))ds (66a) 
Jo 

^v* o(^j^ exp {g{Vf{l - s)) ® exp {g{V)s) ds^ u{dvg{5V)) (66b) 
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hence formally, 

dvgiSV) ^v* o(^j^ cyip{g{Vf{l - s)) ® exp(g(y)s) v{5V) (67a) 

dvg*{A) = cxp(-.g(y)^(l - s)) ® exp(-g(y)s) ds^ v{A) (67b) 

where we have defined the linear operator v : C^^^ — ?> to be the operator that takes a square matrix 
and stacks the columns to create a vector. If we let C^^^ and both be real Hilbert spaces with 
respective inner products {A, B) — '^Tv{A'^ B) and (x, y) — 5R(x'fy), then we can find the adjoint of v: 

2^" 2" 

{v*{^),A) = (x,iy(A)) = ^Y.^kVkiA) = 3? ^ i,+o--i)*2";^»+o-^i)*2"(^) (68a) 

k—1 i^j — l 

= 3?^z.-\x)tA„- =3f^Tr(^.-l(x)^A) = {v-\^),A) (68b) 

for a given x e C^^" and all A £ C^"^^". Hence ~ v^^ is the "matrixization" operator that cuts the 
vector up and arranges the pieces as the columns of matrix. 

By construction, g{V) will be skew-Hermitian and have imaginary eigenvalues in the (imaginary) interval 
(— ITT, itt]. It may be shown that as long as the spectrum of V does not contain the eigenvalue —1 (which 
holds with probability 1), the inverse operator in (j67bl) is well-defined. 

Applying these differentials to our problem, we find that 

grad Jg(J/) = duh* o d^u^g* o grad /(.g(/i(C/))) d^h* o dh{u)9* {g{h{U))) . (69) 

Now, 

dvg* {g{V)) = ^* ° exp ( - g{Vf{l - s)) ® exp ( - g{V)s) ds~^ ,.{g{V)) (70) 

is the matrix solution A to the problem exp ( — g{V)s) A exp ( — g{V){l — s)) ds — g{V). If we assume 
temporarily that A commutes with g{V), then we find that A — Vg{V) is one solution to this equation. Under 
the assumption that the spectrum of Q docs not contain ~1, the inverse operator in ()70p is well-defined, so 
that ([TOl) has a unique solution, and therefore 

dv9*i9{V)) =Vg{V). (71) 

Hence, 

giad Jg{U) = duh*{h{U)g{h{U))) = -Uh{U)g{h{U))W^U (72a) 

= -W\og{U^W)W^U = -U{U^Wlog{U^W)W^U) (72b) 

= -UlogiU'^W). (72c) 
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